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The Suspensions are among the important components of vehicles, 
providing comfort for passengers and protecting the chassis and freight. 
They are normally provided with dampers that mitigate these harmful and 
uncomfortable vibrations. This work aims at investigating the potential of 
the use of the Magnetorheological Damper as a semi-active suspension for 
passenger vehicles to improve their ride dynamics and to, ensure their 
manipulability, and reduce unwanted vibrations levels. To achieve desired 
performances, a hybrid controller based on a backstepping-sliding mode 
control strategy is derived to increase the dynamic performance of the 
automotive suspension component and providing comfortable passenger 
traveling. Numerical simulation is used to test the efficiency of the proposed 
controller. The obtained results show that the proposed backstepping sliding 
mode controller is more efficient and robust against road profile excitation 
and external disturbances in comparison to a classical controller based on 
backstepping and sliding mode control. 
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1. INTRODUCTION 


Nowadays, quality becomes very necessary for people around the world. Especially, comfort in 
vehicles, personal or public, is a priority concern for passengers therefore it is desirable to have a high- 
performance vehicle suspension system. Many vehicle components, such as suspension strut type, rolling 
center location and center of body weight, and steering system, might be potential elements for improving the 
performance of the suspension. Among these choices, suspension absorber studies are several and most 
active because improving performance of damper is the efficient choice to separate undesired vibrations 
generated due road profiles [1]. In the world of moving automobiles, passengers are not fully comfortable 


due to body vehicle vibration. 


As a result, effective regulation of vehicle suspension vibrations is critical for improving driving 
comfort. In general, vehicle suspension is divided into three categories: active, semi-active and passive 
suspension [2]. Passive suspensions, which are often employed in traditional cars, are unable to increase 
driving comfort since their settings are fixed, and they may only perform well on specific road surfaces. 
Active suspension systems can improve riding comfort, but they always demand additional power. 
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Semi-active suspensions systems, on the other hand, feature changing damping coefficients or 
stiffness characteristics and do not require additional power sources. These suspension types are often 
defined by high-level technical professionals and also because of their relatively low frequency, changeable 
stiffness characteristics, and decreased suspension noise [3]. Semi-active suspensions with magneto- 
rheological (MR) dampers have emerged as the finest technical option for improving car suspensions just on 
road or the ride comfort. This seems to be due to benefits like the existence of a compromise in between 
attempting to control requirements for suspension damping for trying to drive comfort and ability to handle 
safety, the suitable control requirement for hysteresis as well as force saturation nonlinearities of MR 
dampers and highway vehicles uncertainty [4]. It has more advantages for obtaining desired performances 
compared with other suspension systems. It has, i.e., semi-active suspension, more advantages for obtaining 
desired performances compared with other suspension systems. For this reason, several researchers and 
academicals have taken more attention to study and develop suitable techniques and methods for semi-active 
vibrations attenuation and suppression. In recent years, experts and researchers have put forward many 
control strategies for vehicle suspension systems, such as PID control, optimal control, fuzzy control, neural 
network and adaptive control, sliding mode control [5]-[9]. 

Prabu et al. [6] have presented a pneumatic actuator control for vehicle suspension system using a 
classical PID controller. The pneumatic actuator model is modeled as the first-order system, and also the PID 
controller design was done using the Ziegler Nicolas tuning approach and optimum control, with comparison 
of the obtained results with passive suspension. Bao et al. [7] created an adaptive SMC controller based on 
fuzzy logic for the active suspension using air spring system. The suggested controller is intended in 
improving active suspension performance by using fuzzy control for dealing with the nonlinearity of system. 
For semi-active air suspension control, an adaptive technique based on neural networks and neural identifiers 
has been developed [8]. The strategy proposed in [8] used both neural networks and identifiers to compensate 
for the time-variation parameters and the nonlinearity of the air suspension system. Chang and Roschke [9] 
devised an alternate model of a MR damper in regards of a neural network. 

A neural network approach with six neuronal inputs or one output to replicate the behavior of the 
MR dampers using the phenomenological model produced data. Wu and Griffin [10] created a relay on-off 
semi-active control strategy with an electro-rheological damper to reduce end-stop effects in suspension 
seats. Choi et al. [11] created a seat suspension system with such a skyhook control method and MR 
dampers. Suitable control techniques for MR systems, as well as design factors [12], have even been actively 
explored to improve the comfort performances with the accuracy and precision under extreme road 
excitations. The backstepping technique is widely used in nonlinear systems control and in suspension 
systems applications [13]-[19]. 

It can ensure worldwide tracking for such nonlinear systems class that can be transformed into the 
parameterized feedback form. Moreover, among novel strategies, the classic nonlinear controller using the 
SMC represents a strong tool for dealing with matched uncertainties and therefore is commonly employed for 
the nonlinear suspension systems [13], [17], [20]-[22]. The sliding mode controller's major shortcoming is 
the undesirable chattering phenomena, which can stimulate high-frequency response thus contribute to the 
instability. Many techniques, such as the combining of backstepping with sliding mode control, have been 
offered to avoid the issue of chattering. These nonlinear controls are intended to increase vehicle agility and 
the driving comfort [13]-[22]. The goal of this research is to create and investigate the vibration signals of a 
MR suspension under various road circumstances employing three nonlinear control systems: backstepping 
strategy, slide mode control, and back-stepping-sliding controlling. The remainder of the study is arranged as 
follows: section 2 develops a MR damper behavior and model. Section 3 includes a quarter car model for 
said semi-active suspension system. Section 4 develops several control techniques for the semi-active 
suspension control based upon backstepping, SMC, and backstepping-SMC design to negate the influence of 
unknown road roughness on the vehicles. The obtained simulation results of these controllers and 
comparisons are discussed in the section 5. Finally, the results are drawn in section 6. 


2. MATHEMATICMODEL OF THE MAGNETO-RHEOLOGICAL DAMPER 

Today, MR dampers have been used for several times to reduce the vibration picks in dynamic 
systems including active suspension control. In several practices cases, it has been verified that undesired 
vibrations of systems can be reduced using MR dampers with known theoretical control strategies [2], [23]- 
[26]. This semi-active device behaves similarly to a hydraulic damper, with the exception that its fluid can 
modify its yield point when such an external magnetic field. Ultimately, MR dampers are strong devices that 
combine the benefits of passive components with the advantages of active control [24], [27], [28]. 

Figure 1 depicts the MR damper arrangement. The damping device is essentially made up of a 
hydraulic cylinder containing magnetizable micro particles flowing in oil-like fluids. A tiny opening on one 
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side of the piston allows fluid to pass through. Once the application of magnetic field is done, the particles 
become polarized and, in a few milliseconds, transition to the semi-solid phase with a high flow resistance. 
The changing magnetic field causes the viscous effects to become adjustable, and as a result, damping forces 
may be created to dampen the existing oscillations [24]-[26]. 


Diaphragm z 
Accumulator Annular Orifice Bearing & Seat 


Wires to 
Coil MR fluid Electromagnet 


Figure 1. The basic configuration of the MR damper 


The mathematic equations of the damper can be used to describe the applied voltage or current or to 
the MR circuit with the produced damping force. Several scholars have constructed statistical models to 
describe the hysteretic behavior of MR damper in the research. Figure 2 depicts the well-known enhanced 
Bouc—Wen model, which describes the characteristics of MR dampers [29]-[31]. Spencer et al. [29] created a 
phenomenological framework wherein the force of the MR damper is described as a function of the spring 
stiffness, accumulation stiffness, viscous damping at a low velocity, and damping or stiffness at high speeds. 
The following is an expression for the forces generated by the MR dampers using the generalized Bouc-Wen 
model: [12], [14], [29]-[31]. 


fur = CY + ky (x — xo) (1) 
Oy = colà — ý) + ko(x — y) + az 2) 
ż = —y|x — ylzlz|""* — BE - ylzl” + AG - y) (3) 
Ý = elas + cot + ko(x + y)] (4) 


Where x, andfypindicate the damper's displacement, speed, supply voltage, and produced force, respectfully, 
z is an evolving variable that whih characterizes the hysteresis behavior of the MR damper force versus the 
damper's movement and velocity. The hysteresis loop's form and size are described by the parameters, and A. 
The device's extra parameters may be represented as a function of applied voltage as follows [29]-[31]: 


A= Ag + Apu; Cy = Cia + Cip ; Co = Coa + Cop (5) 
where u is a variable that defines the dynamic of the system given by: 
u=-t(u-v) (6) 


where v is the control voltage supply at terminals of the driver and Tẹ denotes the first-order filter time 
constant. 

Parameters of the mangnetorheological damper model adopted are given as: 

Coa = 7.84N - s/cm ; Cop = 18.03N-s/om; ky = 36.1N/cm; kı = 8.4N/cm; Cia = 146.49N - s/cm ; 
Cip = 34.622N -s/cm; Qa = 124.41N/cm ; Ap = 384.3N/cm ; y = 136.32cm”? ; 
B = 205.9020cm”? ;A = 251 ;n = 2 ; n = 190571; xo = 09.08cm. 

A sinusoidal signals displacement with the amplitude of 1.25cm and 1.5 Hz frequency is performed 
for varied supplied input voltages to validate the behavior and hysteresis features of the proposed MR damper 
type (0, 0.5, 0.75, 1, and 1.25 V). Figure 3 depict the produced force's temporal response, the pressure against 
MR damping displacement, as well as force-velocity hysteresis. The various plotted charts of the produced 
force clearly show the hysteresis behavior of the MR damper in between time Figure 3(a) and velocity x and 
the damper force fygas shown in Figure 3(b), and the non-linearity characteristics of fyg according to x 
variation is remarkable (see Figure 3(c)). 
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y x 


Bouc-Wen 


Figure 2. The schematic representation of the magnetorheological damper according to 
the Bouc-Wen model 
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Displacement [cm] 
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Figure 3. Produced force of MR damper model on the (a) MR force-time, (b) MR force-velocity and 
(c) MR force-displacement 


3. DYNAMIC MODEL OF MR QUARTER-VEHICLE SUSPENSION SYSTEM 

Despite the fact that the genuine suspension may be treated in a way more sophistical, the suggested 
model simplifies the complicated geometry of the suspension structure [32]. To investigate the robustness of 
the suspension control strategy we suppose a quarter vehicle system as shown in Figure 4 [13], [20], [30]. It 
is a basic structural model since it consists primarily of two subsystems: the wheel mass which represents 
tyre subsystem and also the body mass which represents the subsystem of suspension. These two subsystems 
are connected between them via elastic and dissipative elements and are subjected to external disturbances, 
which are represented by the road profile. From the mechanical model depicted in the Figure 4, the 
mathematic equations of the 1/4-vehicle model are written as follows [13], [20], [30]: 


Msžs = ks(Zu — Zs) + Cs(żu — żs) — Emr 


(7) 


Mužu = hep — Zu) — ks(Zu — Zs) — Cs (Žu — 2) — Enr 
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with z,,z,, and z,denote the vertical displacement of the sprung mass, the unsprung mass and the external 
disturbance due to the road profile, respectively.F,,, the force produced by the MR damper; m,and 
m,,respectively are the car body and the wheel assembly masses.k,and k,, represent the stiffness of the 
suspension and the stiffness of tire.C, is the total damping constant of the suspension. We define the state 
variablesx, the input force vector U and the road profile disturbance Was follows: 


x=[Z, 2, Zu Zy]? =[%1 X2 X3 Xl’, 


then, the 1/4-vehicle model can be expressed in state-space as (8). 


x= Ax + Bu + DW (8) 
Where: 

0 1 
ks G 
Ms Ms 

a 0 0 
ks Cs 
My My 


Figure 4. Quarter-vehicle model with MR damper 


4. CONTROL ALGORITHMS OF VEHICLE SUSPENSION DYNAMICS 

This section proposes semi-active controls of the vehicle suspension dynamics depending on 
multiple nonlinear control schemes. There are three forms of nonlinear control formed: backstepping control, 
slide mode control, as well as a hybridation backstepping-sliding mode control scheme. All of these kinds of 
controls are dependent on the Lyapunov stability criteria, which will be shown in the next subsections. 


4.1. Vehicle suspension control using backstepping strategy 
Backstepping seems to be a systematic and recursive nonlinear feedback control design approach 
ideal for the wide category of feedback nonlinear systems. This technique ensures global pursuit for such 
nonlinear systems category that can be transformed into parametric-strict survey form [13], [14], [33]-[35]. 
The control law's backstepping technique relies on a recursive selection of some relevant functions using 
system states as virtual-control inputs for the under dimension subsystems. Each step of the backstepping 
startegy generates a virtual-control signal, which is expressed in relation of the previous design phases 
virtual-control designs. When technique steps are completed, a feedback disposition for the actual control 
input signal is produced, which meets the desired design target by summing the Lyapunov activities 
connected with each individual design phase [13], [14],[33]-[35]. 
- Step 1: If we opt to reduce the suspension travel, the oscillating subsystem problem will remain 
unaddressed. Because we know that proper controlled variables can help to avoid unpredictable zero 
dynamic features, the regulated output is set to be [36]: 


e1 = X1 —X3 (9) 
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where x; is the displacement of car body, X3 is a filtered version of the wheel displacement x3 : 


Æ x (10) 
where € is a positive constant. 


ĉi =X, - Ž3 
= x3 — €(x3 — 3) 
= x + €(x, — X3) — £61 (11) 


Because of that the error ¢; must be equal to zero, state variable x) is taken as the virtual control to 


make the first subsystem stable. In order to compute the stabilizing function, we defined the following 
Lyapunov function as: 


v =e? (12) 


We differentiate to get: 


V=e,€; 
= €1(X2 + €(x1 — X3) — £e1) 
= —K,e? + (Kye, + x3 + €(x, — x3) — €e;) (13) 


where, K4 is constant gain. The error is achieved if we define stabilization function as: 
a, = —K,e, — E(x — x3) + €e, (14) 


with a, = x; Sand K, > 0 a positive design constant. 
By substituting (14) and (11) into (13) yields 


V, = —K,e? (15) 


- Step 2: now, in the next step we try to make the virtual control xT as desired. So, we define the 


following error signal 
ez = X2 — My, (16) 
and the error equation is: 
ê = —(K, + €)e, +e (17) 
ê, =X. — ay 
= a [k(x = x3) — Cs(%¥2 — x4) — Fa — (Kı (— Ky + £)e, + e2) — €(X2 — X4)] (18) 
To determine the stabilizing function, the augmented Lyapunov function is defined as: 
V =e? +e} (19) 
So, its derivative is defined as follow: 


V = e,€; + e262 = e (—(Ky + €)e1) 
+e, (— [k,(x1 — X3) — C,(X2 — x4) — Fy — (Kı (—Kı + Je, + e2) — €(X2 — tD) (20) 


The control tracking is achieved if one defines the following stabilizing function: 
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Fa = ms(—(K2 + Ky Jen + (K? — 1 + Kye)e1) + ks (%1 — x3) + cs (%2 — x4) (21) 


where F, denotes the desired active force which must be produced by the MR damper (i.e. Fa = Enr). 
Kņis a positive constant, by substituting (21) into (20) yields. 


V, = —(K, + e)e? — Ke3 (22) 


Thus, the system is asymptotically stable. The system is simulated in the Simulink as shown in Figure 5. 


Figure 5. Simulation of a quarter vehicle suspension model using MR damper with backstepping control 


4.2. Vehicle suspension control using sliding mode control 

Many research studies have used the variable structure with the sliding mode, or sliding mode 
control, in the previous decades for non - linear control processes [7], [13], [17], [20]-[22], [37]-[39]. The 
sliding mode controller (SMC) seems to be a non-linear control mechanism that modifies the dynamics of a 
non-linear system, causing it to slide. The SMC is widely renowned for its resilience, high stability, and 
straightforward implementation [40]-[42]. The characteristic of the sliding mode relies on directing a 
system's state trajectory toward this sliding surface and directing it to switch across it to the equilibrium point 
using proper switching logic, resulting in the sliding phenomena. The benefits of SMC control are numerous 
and considerable, including improved accuracy, a very rapid reaction time, and, most importantly, 
resilience [43], [44]. In general, the sliding mode control design procedure can divide in two steps. The first 
step of SMC design involves the selection of the sliding surface (sometimes referred to the switching surface 
or switching manifold S(x) = 0) in the state space where the control undergoes discontinuities [13], [20], 
[21]. This sliding surface is chosen to simulate the required closed loop performances. 

Insecond phase, the control rule should be devised such that such system's state trajectory are driven 
toward and remain on the sliding surface. The reaching phase refers to the system state trajectory before to 
reaching the sliding surface. The generic expression of Slotine and Li's sliding surface is calculated by: 


S(x) = (=< + a)" e(x) (23) 


where e(x) = xg — xdenotes the tracking error, A is a positive constant, xis the variable state. 
The sliding-mode control law u is consisting of two components. One is the equivalent control teg 
and the other is switching control input us, i.e. [7], [20]. 


U = Ueg + Us (24) 
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where Ugg is called equivalent control which dictates the motion of the state trajectory along the sliding 
surface and A is known as switching control part which is given in the basic form as: 


u,=k- sgn(S(x)) (25) 


where kis a positive constant. 
The control input should be designed such that the following sliding condition is satisfied: 


S(x) -S(x) <0 (26) 


using a sign function often causes chattering in practice. Among the solution of this problem is to introduce a 
boundary layer around the switching surface, and the switching component is become as: 


us = k- sa (S09) 27 


where sat(x) is the saturation function defined as (28). 


a peo <é 


sgn (5)/,) S| > ¢ (28) 


wa (8/,) = 


Where € denotes the boundary layer thickness around the sliding surface. The conventional SMC uses a 
control law with large control gains and can produce some drawbacks associated with large control 
chattering. The control-input chattering is undesired and it may excite high-frequency response and lead to 
instability. 


4.2.1. Sliding mode speed controller design 

The goal of the control technique is to create a nonlinear control system based on sliding - mode that 
achieves high accuracy and precision tracking under uncertainties and disturbances. The following stages are 
included in the design of the proposed sliding - mode for the semi-active suspensions control. The following 
is the definition of error tracking: 


e = (xı — X3) 
Then, the general sliding surface is given by: 
d A z 
S= 7% — %3) +1: (1 >) (29) 


a sliding condition must be established in order for the state to travel toward and achieve the sliding surface. 
The suggested sliding condition is [20], [21]: 


S:S<-n-|S| (30) 


where nis a positive constant. 
The damping force is calculated by deriving (29) since it forces trajectory to point toward this 
sliding surface. The systems trajectories, in particularly, adhere to the surface then on surface, i.e. 


S=0 (31) 


From (28) and (29): 
$=ë+1ė 
= ksx; + (Cs = AMsMus/ (Ms + Mus))x2 
(Ks + kusms/ (ms + Mus))x3 + (Cs z AMsMus/ (Ms + Mus))x4 
+ kusXrMs/ (ms ag Mus) E, Fy (32) 
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The best approximation for the desired damping forceF,, a control law should achieve $ = 0 thus; 
ae == kx, + (Cs —Am,m,;/(m, + Mys))X2 
—(k, + kysms/(ms T Mus))x3 T (C; m AMsMus/ (Ms T Mus)) X4 
+ kusxrms/ (ms + Mus) (33) 


In order to satisfy sliding condition, a term discontinuous across the surface is added to the equivalent control 


F;™™ then the desired control force Fa can be expressed as: 
Fy = Fa™ +k sx(5/¢) (34) 


with k > 0. The system is simulated in the Simulink as shown in Figure 6. 
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Figure 6. Simulation of a quarter vehicle suspension model using MR damper with sliding mode control 


Tt 


4.3. The hybrid controller backstepping-sliding mode design 

The control goal is to construct a backstepping sliding-mode control [45]-[49] system to monitor the 
reference trajectory for the outputs of the vehicle suspension illustrated in (9). The suggested backstepping 
sliding-mode control system is intended to meet the goal of automobile body displacement tracking. This 
controller's output is really the reference force which must be delivered by the MR damper, and it is 
discussed in: 
- Step 1: The tracking error is chosen as shown in section 4.1: 


e1 = X1 —X3 (35) 


the derivative of f! is computed as: 
é =% — Xs 
= x3 — £E. (X3 — X3) 
= X + E. (X1 — X3) — E. €61 (36) 
the first Lyapunov function is chosen as: 


V, =e? (37) 


define the first virtual control variable, 
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Xref = P1 = —K1. e2 — E. (X1 — X3) (38) 


then, the derivative of Vi is: 


V, = 210) = (Xp — E ` X3 +E: X1 — E: e) 
= —~(K, +£).e,2 < 0 (39) 


with Kyis a positive design gain, 
- Step 2: We define a sliding surface 


S = xX — Pı (40) 
So, it is easy to find the following expression: 

é, =—-(K, + £)e +S (41) 
The second Lyapunov function is chosen as: 


V =V, +25? (42) 


where S is the sliding surface. So, the derivative of V, is expressed as follow 
V,=V,+5-S (43) 


To ensure that VÝ, < 0 we take $ = —Kş sgn(S) — K - S ; with Ky and K are positive gains. 
Differentiating the (40) which yields. 


S=%,— Gy = (ks (3 = x1) + Coq — x2) + Fa) — Gr (44) 
Finally, the control law resulting by sliding mode is of the form: 

k= m,(—Ky sgn(S) — K -S) — ks(x3 — x1) — Cs(x4 — x2) + MG, (45) 
then, 

V, <0 (46) 


therefore, e} and e,will converge to zero. The stability of the proposed backstepping sliding mode control 
system can be guaranteed. The system is simulated in the Simulink as shown in Figure 7. 
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Figure 7. Simulation of a quarter vehicle suspension model using MR damper with hybrid controller 
backstepping-sliding mode 
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4.4. Semi-active control system 

The semi-active control system consists of two controllers. The first one is the system controller 
which produces the appropriate active damping force (F,) corresponding to the suspension responses due to 
the road profile excitation, while the role of the second controller is the adaptation of the input voltage level 
which keeps the MR damper force on the desired active damping force (F,). In this paper, the semi-active 
controller is essentially composed of backstepping controller, sliding mode control or backstepping sliding 
mode control (system controller) and a function unit which used for converting the desired control force to an 
input control voltage (damper controller). The usually adopted algorithm for the conversion of the force to a 
adjusted voltage for the conversion is a clipping algorithm [29], [32]: 


v= v{(Fa E Frid) i Fnrłmax (47) 


where Vmax is the maximum applied voltage, H{-} is a Heaviside step function,Fyp is a measured MR 
damper force, and F,,,denotes the desired active control force produced by the system controller 
(backstepping control/SMC or Backstepping-sliding mode control (B-SMC)). 

The structure of the closed-loop semi-active controller using the proposed nonlinear controllers for 
vehicle suspension with an MR damper is depicted in Figure 8. In the first step, the actions specified control 
force is computed using one of the nonlinear controllers based on vehicle suspension system outputs. 
Because the control input v directly represents the damping forcefyp, the Heaviside step function is 
employed to adjust the input voltage control according to the desired damping function and the real force 
realized by the MR damper. 


Control force 


Figure 8. The block diagram of a semi-active control system for vehicle suspension with MR damper uses the 
proposed nonlinear controllers 


5. SIMULATION RESULTS 
The quarter vehicles suspension structure is subjected to two types of road stimulation in order to 
assess the performance of the semi-active control with MR damper utilising the suggested nonlinear 
controllers. The bump in Figure 9(a) is the very first stimulation signal studied, and it is typically used to 
highlight dynamic response properties. The bump excitation is described by: 
A = A,,(1 — cos(wt))1 <t <2 (48) 


Zy = Ootherwise 


where w = 2mfdenotes the frequency [rad/sec], A,,represents the bump height. The second road excitation 
tested is a step excitation in Figure 9(b) with an amplitude of 0.2 m. These types of road excitations, 
illustrated in Figures 9(a) and 9(b) are applied to the vehicle suspension system as shown in Figure 8. 

The following are the different parameters of the suggested adaptive sliding mode control that are 
regarded to produce high transient control stability and performance which are used for control loop: 
K, = 450, K, = 15, k = 12.5, § = 0.05, K; = 7.5, K = 2.5. The simulation parameter values of the vehicle 
suspension system are listed in Table 1. 
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Bump input [m] Step input [m] 
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Figure 9. Type of road excitations for simulation: (a) bump excitation and (b) step excitation 


Table 1. Parameter numerical value of the quarter vehicle model 


Symbol Value Description 
Ms 315 kg Sprung mass 
Mus 37.5 kg Unsprung mass 
K; 29500N/m Suspension linearised stiffness 
C; 15 00N/m/s Suspension linearised damping 
Kus 210000N/m Tyrestiffeness 


Figures 10(a)-(c) show the vertical displacement of an automobile body, acceleration of the sprung 
mass, and vertical displacement of the unsprung mass utilizing nonlinear controllers (backstepping control, 
sliding mode control, and backstepping-sliding control) for the bump road excitation. In contrast to the 
passive suspension, the outcomes shown in Figure 10 clearly demonstrate that using the suggested semi- 
active suspension system can obtain excellent performances in which the undesired vibrations due to road 
excitation have really been clearly diminished (lower amplitude for the vehicle acceleration, vehicle 
displacement, and wheel displacement). It definitely demonstrates the efficacy and durability of the 
recommended semi-active vehicle suspension control with MR dampers for induced vibrations mitigation. 

Figures 11 (a)-(c) depicts the acquired results for vertically vehicle displacement, vertical 
accelerating of sprung mass, and wheels vertical displacement for step road excitation. According 
to simulation results, the systems displacement and accelerating responses are significantly decreased, and 
the steady-state duration of controller system responds is significantly quicker than the passive system. We 
see that the displacements as well as accelerations have been reduced significantly, confirming that the 
performance have been enhanced. Thus, it can be stated that high passenger comfort may be gained by 
minimizing body movement as much as feasible, as demonstrated by the nonlinear semi-active control. Based 
on the dynamic response characteristics, we can conclude that the suggested semi-active control with the MR 
damper based on backstepping-sliding mode approach outperforms the other investigated nonlinear controls 
(i.e. Sliding mode control and backstepping control). 

Figure 12 depicts the applied voltage level to the semi-active suspension's MR damper. We can see 
that the Voltage level is constantly switching between zero and the maximum voltage based on differences 
between produced forces and clipped-control. Figure 13 depicts a contrast of the real damping force, desired 
damping force, and active force. Figure 13 means that the actual damping forces generated by the MR 
dampers can following the reference damping force with precision and decrease tracking error, demonstrating 
the usefulness of the clipper optimum control utilised in damping force control. 
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Figure 10. Simulation findings for semi-active vehicle suspension control with MR damper under bump road 
stimulation: (a) vehicle body displacement responses (b) vehicle body accelerating responses, and (c) wheel 
displacement responses 
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Figure 11. Simulation findings for semi-active vehicle suspension control with MR damper under step road 
stimulation: (a) vehicles body displacement reactions, (b) vehicle body accelerating responses, and (c) wheel 
displacement reactions 
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Figure 12. depicts the input voltage provided to the Figure 13. A comparison of the real damp force and 
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6. CONCLUSION 

To reduce undesirable vibrations caused by road irregularity, a semi-active control based on a 
clipper optimum control method for MR damper is suggested in this study for a vehicle suspension system. 
To improve the performance of the vehicle MR damper suspension system, the suggested control strategy 
includes the creation of a classic control scheme, a back stepping control approach, and a hybrid back 
stepping-sliding phase controller. The nonlinear control techniques and the clipper control algorithm are used 
to generate the active damping power and the MR damper driver's modified input voltage. In the first stage, a 
generalized model of the MR damper based on modified Bouc-Wen model is provided. Next, for a semi- 
active controlling, 3 types of nonlinear control methods, backstepping, slide mode, and backstepping-sliding 
authority, are created and used to the quarter-vehicle rear suspension. 

Numerical methods utilizing MATLAB/Simulink programs are used to evaluate the effectiveness of 
the suggested control techniques for various forms of road excitations. It has been quantitatively 
demonstrated that the semi-active controller performs well and can significantly minimize or suppress the 
produced vibrations. Moreover, we observe that the hybrid backstepping-sliding mode controlling is more 
durable and effective than the other nonlinear systems. Thus, the obtained findings indicate that proposed 
control mechanisms connected with the chosen MR suspension system may greatly increase passenger 
comfort. 
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